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Effective-one-body Hamiltonian with next-to-leading order spin-spin coupling 


Simone BalmellQ and Philippe Jetzer 

Physik-Institut, Universitat Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland 

We propose a way of including the next-to-leading (NLO) order spin-spin coupling into an effective- 
one-body (EOB) Hamiltonian. This work extends [S. Balmelli and P. Jetzer, Phys. Rev. D 87, 
124036 (2013)], which is restricted to the case of equatorial orbits and aligned spins, to general 
orbits with arbitrary spin orientations. This is done applying appropriate canonical phase-space 
transformations to the NLO spin-spin Hamiltonian in Arnowitt-Deser-Misner (ADM) coordinates, 
and systematically adding “effective” quantities at NLO to all spin-squared terms appearing in the 
EOB Hamiltonian. As required by consistency, the introduced quantities reduce to zero in the test- 
mass limit. We expose the result both in a general gauge and in a gauge-fixed form. The last is 
chosen such as to minimize the number of new coefficients that have to be inserted into the effective 
spin squared. As a result, the 25 parameters that describe the ADM NLO spin-spin dynamics get 
condensed into only 12 EOB terms. 

PACS numbers: 04.25.-g, 04.25.dg 


I. INTRODUCTION 

Thanks to the LIGO/Virgo network of second- 
generation ground based interferometers, a first direct 
detection of gravitational waves (GW) is expected to oc¬ 
cur in few years [T]. Furthermore, in the next decades, 
space-born GW detectors [2] (such as the planned eLISA) 
will open an entire new window to astrophysics, and allow 
tests of General Relativity at an unprecedented level [3]- 
Both types of detectors rely on coalescing (and possible 
spinning) black hole binaries as a primary GW source. 

In order to extract the GW signal from the noise, very 
precise waveform templates need to be constructed. Gur- 
rently, the most accurate description of coalescing black 
holes is provided by numerical relativity (NR) (see e.g. 
[1H7] for some recent advances). However, the full param¬ 
eter space is too large (especially in the case of nonzero 
spins) for being densely covered by NR simulations . 
This is the main reason why semi-analytical methods, 
by far less computationally expensive, can turn out to be 
very useful. At the present time, the effective-one-body 
(EOB) approach (we refer to [S] for a general review) is 
the only semi-analytical method that has been able to ac¬ 
curately describe the complete waveform of a coalescing 
process. After years of constant development [5H21] the 
EOB has now reached an excellent agreement with NR 
waveforms in the case of nonspinning binaries [ZlEl]. A 
lot of effort has also been put into the modeling of spins 
[23H2ni : leading to a good overlap with NR waveforms in 
the case of nonprecessing (aligned or anti-aligned) spins 
[II ISHl |40]. By contrast, the description of processing 
spins still needs some work before reaching a comparable 
performance |39j , and is currently one of the most urgent 
tasks. In view of this, it may be crucial to incorporate 
more analytical information from the spin-orbit and spin- 
spin Hamiltonian computed within the post-newtonian 
(PN) theory. 
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In Ref. [35], a possible way of including the next-to- 
leading (NLO) spin-spin coupling HT] |32] (see also [331 - 
l48] ) into the spinning EOB model of Refs. [^[27l[34] has 
been exposed for the special case of equatorial orbits and 
nonprecessing spins. A general inclusion of NLO spin- 
spin effects for fully processing orbits would be a neces¬ 
sary step for extending the reliabilty of EOB waveforms 
to a significantly larger portion of the parameter space. 
The present paper aims at filling this gap, providing a 
possible implementation of the missing terms. 

Recently, an improved (and calibrated) EOB model 
for spinning binaries has been proposed iOl, where the 
NLO spin-spin coupling is only incorporated for the case 
of circular orbits. The present paper could be a first step 
for developing, at a next stage, a more general version of 
that improved model. 

The paper has the following structure: in Sec. 0 
we summarize the main concepts and the formalism of 
Ref. [36]. In Sec. (HI|, which is the central part of the 
paper, we propose a way of including the NLO spin-spin 
terms and show the explicit result, both in a general form 
(with a given number of free gauge parameters) and in 
a gauge-fixed formulation. Finally, Sec. (IV) discusses a 
second, slightly different way of incorporating the wished 
spin-spin terms. Both approaches are compared plotting 
the (gauge invariant) angular frequency and binding en¬ 
ergy at the last stable orbit (LSO). The plot also shows 
the prediction of the calibrated models of Refs. [351135] . 

Throughout the paper, we use geometric units with 
G = c = 1. When writing formulae in a PN ex¬ 
panded form, however, we will reintroduce the usage of 
c, prepending a factor (1/c) with the mere purpose of 
labeling the PN order n. 


II. SUMMARY OF THE PREVIOUS WORK 

In this section, we outline the method followed in 
Ref. [36], which forms the basis of the current paper. We 
work in EOB coordinates, with R being the radial coor- 
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dinate, n the unit radial vector and P the momentum 
vector. We will often use the rescaled variables r = R/M 
and p = P/fi. Here, M = mi + mi is the central EOB 
mass (mi and m 2 being the individual masses of the two 
black holes), and p = mim-ijM is the reduced mass. 
Moreover, we denote by = p,/M the symmetric mass 
ratio. 

The starting point is the EOB model of Ref. [55] (which 
includes both spin-spin and spin-orbit coupling at lead¬ 
ing order (LO)), together with its extensions to the NLO 
m and to the next-to-next-to leading order (NNLO) [Mj 
spin-orbit coupling. As already mentioned, in Ref. [35] 
the NLO spin-spin Hamiltonian in ADM coordinates has 
been reformulated and inserted into this EOB model for 
the special case of two black holes whose spins are aligned 
(or anti-aligned) with the orbital angular momentum. In 
particular, it has been shown that it is sufScient to re¬ 
place the spin parameter og of the effective metric (see 
for instance Eqs. (4.7) and (4.8) of Ref. |35j), when¬ 
ever it appears as a second power, by a new, effective 
squared spin parameter. Using the dimensionless nota¬ 
tion xo = 0 , 0 /M, this prescription takes the form 

xg^(X^)eff = xg + AXeV (2.1) 

We recall that xo is a combination of the dimensionless 
spin parameters xi and X 2 (Xa = Sa/m^) of the two 
bodies: 


Xo 


mi 


Xi + 


m 2 


M 


X2- 


( 2 . 2 ) 


The additional term Axeg is of fractional IPN order with 
respect to Xo and carries the information for reproducing 
the correct NLO spin-spin coupling. It reads as 


AXeff 


1 / 2 I CiiN 2 

( 2 I O22 \ 2 

+ \ a22P + —j X2 


ai2P 


Cl2 \ 
r / 


X 1 X 2 


(2.3) 


The calculation of the coefficients Uab and Cat is the 
main result of Ref. [35], given by Eq. (5.12) there. All 
of them vanish in the test mass limit v ^ 0, consistently 
with the requirement that the EOB metric must reduce 
to the Kerr one. 

PN results in ADM coordinates can be included into 
an EOB model after suitable canonical transformations. 
We denote here by G™ and G™ the generating func¬ 
tions of the corresponding purely orbital and spin-spin 
transformation, respectively. The procedure for trans¬ 
forming the NLO spin-spin Hamiltonian from ADM into 
EOB coordinates is made of three stepsQ 


^ Here, we treat the PN expansion under the assumption of rapidly 


1) A purely orbital transformation 

J^NLO/ _ £^NL0(ADM) |(;j1PN ^LO(ADM) | 

2) A LO spin-spin transformation 

rrNLO'/ rjNLO/ i //^LO rrlPN/'l 
^ss = ^ss + l^SS 

where 

jjlPN/ _ ^IPN(ADM) |(^1PN ^N(ADM)| 

3) A NLO spin-spin transformation 
jjNLO/ff _ ^NLO// _|_ |giNLO ^N(ADM)| 


(2.4) 


(2.5) 

( 2 . 6 ) 


(2.7) 


Notice that, for consistency, the transformations must 
be performed in a well-defined order. As indicated above, 
we make first use of the orbital transformation Go, and 
then of the spin-spin transformation Gss Notice that, 
in Ref. [34] , the spin-orbit generating function Gso is also 
applied after Go- By contrast, an evaluation order pre¬ 
scription between spin-orbit and spin-spin transforma¬ 
tion would first be necessary when taking into account 
contributions that are cubic in the spins. 

The final Hamiltonian must be equal to the corre¬ 
sponding term arising from the PN expansion of the EOB 
Hamiltonian. Since, for the moment, this is only true in 
the spin-aligned case, we are just allowed to write: 


rrNLO/// _ ^NLO(EOB) 

.^ss,al ~ ^ss,al 


( 2 . 8 ) 


The NLO spin-spin transformation G^lo = G^lO/^ 
is given by 


ANLO _ • P) 

'-^ss,al ^6 j. 
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(2.9) 


rotating black holes {Sa = -^Xa, with |xa| $ 1), which as¬ 
signs a well-defined PN order to the spin-dependent terms. As 
a consequence, throughout this paper, Gss and Hss are of 2PN 
order when labeled with “LO”, of 3PN order when labeled with 
“NLO”, and so on. 

^ Since the set of canonical transformations carries a group struc¬ 
ture, the successive evaluation of Go(r,p^) and of Gss(i*^p^^) 
is a canonical transformation itself (with generating function 
Gss,o(r,p^0 = Go(r, p^) -h Gss(r^p^0 - ■ p^). Despite tak¬ 

ing a unique generating function would avoid the necessity of 
fixing an evaluation order prescription, we prefer here to use two 
separated transformations, so as to mantain the continuity with 
respect to Ref. ESj. As a second reason, the transformation of 
the Hamiltonian is more easily calculated here than for a sin¬ 
gle generating function, since in that last case effects quadratic 
in the transformation must be considered (see e.g. Eq. (6.9) of 
Ref. [9]). 
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Notice that, in view of the generalization to processing 
orbits, we have written the individual dimensionless spins 
Xa — as vectors. The coefficients Uab and ^ab can 

be found in Eq. (5.13) of Ref. [5S]. It is also provided 
an expression for in the test mass limit (assuming 

TOi > TO 2 ): 


limG^LO 

[/—J-O 


1 


\{xl + {n-Xif){n-p) 
{P-Xi){n-Xi) ■ 


( 2 . 10 ) 


The purpose of this paper is that of generalizing the 
prescription (2.1) to the case of general orbits, i.e., when 
the scalar products (n • Xa) and {p ■ Xa) (a = 1) 2) can¬ 
not be set to zero. 


III. INCLUDING NLO SPIN-SPIN TERMS 
INTO THE EOB FOR GENERAL ORBITS 

A. The prescription 


In Ref [36], the EOB metric is written in Boyer- 
Lindquist-like coordinates. When spin precessions must 
be taken into account, however, it is necessary to switch 
to another system of coordinates. Following Ref. [26] 
(and reformulating the angular variable 9 according to 
oo cos{9) = {n ■ ao)), we can write the effective metric in 


Cartesian-like coordinates, 

geff - (^ ■ Qo)^ - ^ ^ (3.1) 

(3.2) 

g% =^(^Ann^n^ + {5^^ - nV) (3.3) 

_^, (aoxn)>oxn)^ \ 

At ) 

where f?' = R? + [n ■ ao)^- and A/j can be found 
e.g. inEq. (4.9) of Ref. |3Bj (notice that they both depend 
on the spin through a term ^ Oq). The effective Hamil¬ 
tonian (that we denote here as “old”, in order to avoid 
confusion with the modified version that is presented in 
this paper) takes the form 

Bfi = yM"+7bp,P^ +Q4(P,), (3.5) 

with a quartic-in-momenta term Qi{Pi) |llLI26j and with 


1 

y-sSS 

(3.6a) 


(3.6b) 

9eS+ ^2 ■ 

(3.6c) 


AHso has been introduced to describe higher-order spin- 
orbit couplings. It is defined in terms of the gyro- 
gravitomagnetic ratios and g'gl, see Eqs. (4.15) and 
(4.16) of Ref. Elj. 

With the reduced quantities Xo = clq/M = {rriiXi + 
m 2 X 2 )/M, H = H/g, A^ = A^/M and A* = At/M, we 
can write, more explicitly: 


AHso + P^Pi 


a 


+ (^(j^gf + g'^l^ inxp)-xi+ (j^^gf + gf^'^ {n x p) ■ X2^ 


(3.7a) 


(3.7b) 




P 


-1 I (n.p)^- 1 (2r^ 


-^t + xl + {n- Xo)^ ) {{n X p) ■ Xo) 


, (3.7c) 


where 

+ Xo)^ - (Xo “ (« ' Xo)^) ■ 


The spin-squared term {{n x p) ■ xof', generated by the 
contraction of Y^PiPj > can be expressed through the sim¬ 
ple scalars p^, (n-p), Xg, (n • xo) and (p • Xo) according 


(3.8) 



















4 


to 


{{n X p) xof = {p^ - {n ■ pf) (xo - (« ' Xo)^) 

- [{P ■ Xo) - {n ■ P){n ■ Xo)f ■ (3.9) 
Let us now do some considerations: 


The prescription (2.1) acts selectively - leaving all 
non-squared spins untouched - and cannot be truly 
considered as a redefinition of the effective spin of 
the EOB metric. It is, rather, a direct modification 
of the effective Hamiltonian. 


One can say the same for the inclusions of spin-orbit 
terms done in Refs. [551 HZl [M]. Adding the spin- 
orbit coupling requires a modification of all terms 
in the metric that are linear in the spin - or, equiv¬ 
alently, a direct modification of the Hamiltonian 
through an additional quantity AHso (Eq. (4.16) 
of Ref. |57]b 


ii) Changing the Hamiltonian itself rather than the 
metric is of course not unreasonable. Since the mo¬ 
tion of a spinning particle is non-geodesic, there is 
no reason to believe that the dynamics of spinning 
bodies can be accurately described by geodesics in 
an effective metric. One should in principle not 
worry about intervening on the structure of the ef¬ 
fective Hamiltonian itself. 


In view of these remarks, and noting, in addition, that 
the effective Hamiltonian depends on the spin squared 
only through the scalars Xoj (’^ ' Xo)^) {p ■ Xo)^_^d 
(*^ ■ Xo)(p • Xo)) we see a natural way to generalize (2.11. 
We treat the spin differently whether it is contracted with 
itself, n or p, and propose the following type of replace¬ 
ments in Eq. (3.71: 


Xo 


(x')eff =Xo + ;:^ 


Jx)^2 I Jx)^2 I Jx)^ ^ 
^11 Xl + 2^22 X2 + H2 Xl X2 


{n-Xof 

{p ■ Xof 


(« • x)eff = (« • Xo)^ + 


(p ■ x)eff = (p ■ Xof + -y 


• Xi)^ + 42 ^(« • X2)^ + z[l\n ■ Xi)(n • X2) 
Zn{p ■ Xi)^ + z!f^{p ■ X2f + z[l\p ■ Xi)(P • X2) 


(n ■ Xo)(p ■ Xo) ((w • X){P ■ X))eff = (»^ • Xo)(p ■ Xo) + ^ 


z[T\'n • Xi)(P • Xl) + ■ X 2 )(p ■ X 2 ) 


+ \ • Xi)(P • X2) + 4 r^(P • Xi)(« • X2)) 


(3.10a) 

(3.10b) 

(3.10c) 


(3.10d) 


where 


— I , . ( 


^ab 


= a. 


VP^ + b^ab\n-py 



(3.11) 


the symbol (x) corresponding to (x), (n), (p) or (np). 
Recall that the functions A,, and At depend on Xo? &iid 
thus need to be transformed according to ( 3.10a |. The 
effective Hamiltonian that results applying (3.10) onto 
H°^ will be simply denoted as Heg. 

The coefficients and are already known, 


I 

their explicit expression being given by Eq. (5.12) of 
Ref. [35] (where the label (x) had not be en used). In 
particular , sin ce the vanish, (3.10a) is consistent 

with Eq. (2.3). 


Determining the z^’s from the Hamiltonian in ADM 
coordinates cannot be done without finding, simultane¬ 
ously, the generating function of the corresponding 

canonical transformation (see Eq. (2.7)). We are looking 


for a sufficiently general ansatz that implements its al¬ 
ready known test-mass limit (2.10), and that mantains, 


in addition, the symmetry under exchange of the labels 1 
and 2. The searched canonical tranformation may have 
the following form: 
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ANLO _ 
— 


(n-p) 


+ {n-p) 


Cn^ - X? + fCn^ “ ^ 1 + '^11^ 


C22^ - ^) X2 + fc22^ “ ^ ■ ^2)^ + ^22^ (P • X2)^ 


+ {n-p) d 2 ^(Xi • X 2 ) + C!i^\n ■ Xi){n ■ X 2 ) + 5^iHp ' Xi){p ' X 2 ) 
+ \ (dr^(« ■Xi){p- X 2 ) + C27\p ■Xi){n- X 2 )) |, 


+ ( dr^ + - ) in-Xi){p-Xi) 


^ ) {n ■ X 2 ){p ■ X 2 ) 


(3.12) 


r 


where 


Aixj _ 
^ab 


«SV + -5;j'(n-pf+ ^ 


W' 

Tab 


(3.13) 


for (x) = (x), (n) and (np). The instead, con¬ 

stant coefficients. Notice that a canonical transformation 
of this type applies an infinitesimal rotation on the spins 
according to 


Xa = Xa 


dG 

dXa 


X Xa 


In our specific case, the spins are left invariant under 
the constraint of aligned spins and equatorial orbits: 


1 dXa 


X Xa 


= 0 . 


al 


formation given by (3.12), which is a consistency require¬ 


ment for the approach we are following. 


B. The general solution 


In order to determine all coefficients, we first explic¬ 


itly calculate the transformations given by Eqs. (2.4)- 
(2.7). All needed expressions are already collected in 
Ref. [35], and specifically: c^n be found in 

Eq. (2.3) there, 7 j 1 PN(adm) ^ 2 . 4 ); in 


Eq. (2.7), H, 


NLO(ADM) 


in Eq. (2.9); in Eq. (3.3), 


and G^s° in Eq. (3.5). We need to rearrange two of 
these formulae, namely ^nlo(adm) G^g°, that had 
been originally written in a form that is not suitable for 
our purpose. In Eq. (2.9b) of Ref. [35], it appears the 
scala r {{p x xi) ■ n) ((p x X 2 ) • n). Using the identity 
(3.9) (but with the vector xi + X 2 instead of Xo) it is 


Thus, nonprecessing orbits are preserved under the trans- easy to show that it can be decomposed as 


([P X Xi) • «)((P X X 2 ) • n) ={p^ - {n-pf){xi ' X 2 )-p^{n ■ Xi)(« ' X 2 ) - {p ■ Xi)(p- X 2 ) 

+ {n-p){{p ■ Xi)(« • X 2 ) + [n ■ Xi)(P • X 2 ))- 

Eq. (2.9b) of Ref. |35| then becomes 

'1 i/' 


TyNLO 

^SiS2 


+ 0 Pdxi • X 2 ) + (1 - 0 (n • pfixi ■ X 2 ) 


+ (^ + 0^ Pd« • Xi){n ■ X 2 ) + 0(n • pf{n ■ xi)(n • X 2 ) 

+ (^ + ^) (P-Xi)(p- X 2 ) - (1 + ^ {n ■ p){n ■ xi)iP ■ X2) 


- ( 1 + ^- («--p)(p-Xi)(w-X2) 


6 (xi • X2) - 12(n • Xi)(« • X2) 


(3.14) 


(3.15) 


Furthermore, using basic vector identities, Eq. (3.5) of 
Ref. |55| is simplifed as follows: 
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1 1 
2 

+ (Xo ■n){rxp)- {xo x n) | 


[Xo - (Xo ■ nf] {r ■ p) 


1 1 ■ 
c'^ 2r . 


{n-p)xl-{n-Xo)ip-Xo) ■ (3.16) 


one), with 72 variables. This means that, if the system 
admits a solution, there will be at least 15 undetermined 
variables, that, as we shall see, will play the role of gauge 
coefficients. Notice that, because of the symmetry under 
exchange of the particle label 1 and 2, the system can be 
reduced to 39 equations and 50 variables. The general 
set of equations is solved by: 


After all transformations (2.4|-(2.7), the resulting 
jjNhOiii be equated to the corresponding term 

^NLO(EOB) ebtained by a PN expansion of the EOB 
Hamiltonian 


Heob 


A 


l + 2v 



(3.17) 


The expansion can be done simply replacing r —>■ r/e^, 
p ep and performing a Taylor series in the small num¬ 
ber e. is then defined as the part propor¬ 

tional to e® which is quadratic in the spins. Finding a 
solution for the equation 

(r,p) = i/NLO(EOB) 


is equivalent to solving an inhomogeneous system of 57 
linear equations (18 for the spin(l)-spin(l) combination, 
18 for the spin(2)-spin(2) and 21 for the spin(l)-spin(2) 


Ax) 


(n) 

ail = 

„,(x) _ 




0^12 — 


lli^^ ^ mi 
32 

0 



0 

0 

lli/2 1/ 


0 

2 

0 

0 . 


TO2 

mi 


(3.19) 


(x) _ _ 
^11 — 

b[f = 0 

Ax) _ 

Cl 1 — 


mi 

16 2 m2 

29i^^ mi 

16 2 m2 


(y) 

=-v- 


lli/2 








Q.11 — 


19v 39v^ f V 


ISZx' \ 77T-2 (n) (np) 


A") 

"11 


bv 15z/^ / bv Ibv'^ 


- 57^?^ - 

mi 


Hi/ 


Cii — 


+ ^ +1-2 


7iA 

4 


m2 

mi 


(") 

■7ii 


(np) 

111 


®12 “ 


L/i o — 


— 

^12 “ 


r, I ^ f I ('^P) 

- 2 !^ -^ + 7 i 2 - 2 (,“21 + “12 


-5iy- 
-2v 


4 
15i^^ 
2 

2 


) 

54?-1 (/Sir*+dr’) 


(n) 

7i2 


>7ir’ 


7lr’) 
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ip) 

«ii 

bip) 

„ip) 


ip) 

«12 

biP) 

^12 

Ip) 

<-12 

,("P) 

^11 

!,("P) 


Anp) 

-11 

,("P) 

•^12 

iXnp) 


Jnp) 

-12 

,(«P) 

-^21 

!,("P) 


„(«P) 

-21 


2(aS7’ + «l?) 

2(/3(7) - 2 SiP ^) 

47)+47)+25(^2^ 

/37p)+77)_45W 

-i2 + i.2+77)+77)_25W 

4/3(7) 

^ ^ + ( 4 . + ^ 1 (-87(?) + Salr^ + 8/3^^^ + 

27(7)-/3(7)) 

4/3(7) 

_^ _ 12^3 - (5122 ^ g^3) ^ _ (2^2 ^ ^ _ 27") + 2a(7) + 2/3(7) + 37^7^ + 25^^) 

2(47)-/37^)) 

4/37") 

_^ _ 12123 _ (5^2 ^ g^3) ^ _ (2^2 + 6123) !!Ji _ 27^) + 2a7") + 2/37") + 377'’^ + 2^4 


( 3 . 20 ) 


As already mentioned, the coefficients 022, /322, 722, 
^22, ^22 and C22 directly follow from the solution above 
exchanging the particle labels 1 and 2 . In regard to this 
point, it is worth discussing the Kerr limit 12 —>■ 0 . In this 
case, indeed, one has to choose which one of the two spins 
vanishes - thus somehow breaking the symmetry between 
them. For 12 —>• 0 (and m2/mi —?► 0 ) all coefhcents a, j 3 
and 7 must vanish, with the only exception of 77^ > 77^ 
and 77^^- The reason of that lies in the form of the 
canonical transformation ( 3 . 12 ) and its limit ( 2 . 10 ). The 


price of having enforced the formal symmetry between 
spin(l) and spin(2), by adding i2-independent terms also 
to the spin(2)-spin(2) contribution of is that of 

generating coefficients 722 that do not tend to zero in 
the Kerr limit (and specifically: 7^) — 1 / 2 , 77 ^ 1/2 

and 77^^ ~lj which can be easily verified taking into 
account that vm\/m2 —> 1). 


By contrast, all coefficients a, b and c vanish for v —>■ 
0 , as required by consistency. In particular, the non¬ 
zero limit of 77\ 77^ ®'' 3<1 77^^ responsible for the 

convergence towards zero of ; <^7^ > ^22 c7^^, 

which might not have seemed immediately obvious from 
Eq. ( 3 . 20 ). Alternatively, to make this convergence more 


47^77^ 


~(n) _ (n) 

722 = 722 


12 

y 


m2 j V 12 
mi \ 2 4 


~(np) _ (np) , 

722 = 722 + 1 > 


( 3 . 21 ) 

which absorb the_j2-independent terms in the spin(2)- 
spin( 2 ) part of ( 3 . 12 ), and satisfy 47 ^ 4754^^ 0 - 


Now we can do the following reformulation: 




2iy 


Goo - 


&(”) - - 
*^22 ~ 


M _ _ 
<-22 ~ 


+ 


5j2 

y 


12 

2 ’ 

isy 

12 

2 ~ 


15j2^ 

m2 

mi 

7i2^\ m 2 

4 J mi 


^2 ~ (n) inp) 

- 1 " 722 “ «22 


- 547 - 47 ) 


-in) 

■722 


~ ("P) 

722 


4) = y _ 2:2^^ 


jnp) 

<-22 


2 

34 

4 


+ 247 ^ 


mi 


247 )- 257 ) 


- 4 i 2 + 


3i2^ 


y - 247 

mi 


247 ) 


347 ^ 


2 S. 


ip) 


( 3 . 22 ) 


explicit, one could redefine 


We stress that, in this case, the formal symmetry with 
respect to the corresponding on , bn and Cn is not di¬ 
rectly visible. 











C. Gauge fixing 


The general solution (3.20) contains 18 gauge param¬ 
eters, namely 7 ^"^ S[f , a[Y\ 7 ^ 7 ^ 722 ^ ^ 22 ^ 

(np) ^(np) (np), , (n) Ap) (np) (np) Anp) 

^22 ’ A ^22 ’ /22 ’ ai-iu. ^-^2 5 ^12 ’ ^12 ’ ^21 ’ A ^12 ’ 

/32i^\ 7 i 2 *^\ 721 ^^- An appropriate choice of them can 
be used to simplify the 30 coefficients and 

(for (x) = (n), (p), (np)), that are not uniquely deter¬ 
mined. One could try to impose the so-called Damour- 
Jaranowski-Schafer (DJS) gauge (see e.g. Refs. [27l IM] 1 
which would consist in eliminating those new terms pro¬ 
portional to p^, i.e., all coefficients We cannot, how¬ 
ever, make all of them vanish: there is no way, within the 
method we have followed, to impose the DJS gauge. No¬ 
tice that this does not mean that the DJS gauge is a bad 
choice in general: in the spin-orbit sector, it is a very 
useful gauge, independently of how the spin-spin sector 
looks like. Moreover, the impossibility of imposing it in 
this paper is strictly related with the type of prescription 
we have followed. It would be interesting to investigate, 
in a future work, if a different method for including NLO 
spin-spin terms would allow to impose a DJS gauge in 
the spin-spin sector too. 

Here, we nevertheless wish to give an example of gauge 
fixing, and choose an alternative approach. We can, for 
example, search a gauge for which the maximal number of 
coefficients can be set to zero. One finds that at most 24 
of them can vanish, the remarkable fact being that this 
happens for a unique gauge choice, and namely when 


(n) (n) (■ 

11 I <^11 j ^22 , 1-22 J 




a 


(n) 

12 


the non-zero coefficients are a 

and Cj 2 ^. In this case, the gauge coefficients are fixed as 
follows: 


(n) V 

= 4 

= 0 


m2 

mi 


Ajip) 

^11 


= 0 


7ii 


= —2v - V — 

4 mi 


=0 

(«) 3i/2 

7l2^ = --^ 


= 0 
= 0 
(np) ^ 

721 =y 
= 0 
= 0 
(np) ^ 

712 =y 

= 0. 


^2^ 

mi 


2 mi 
y - 
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(3.23) 


We may also write 


~ (") 
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~ ("P) 
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V 
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mi 
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mi 


(3.24) 


The coefficients of the effective spin squared take the fol¬ 
lowing, remarkably simple form: 


7x) _ 
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I 8 2 m2 


7;’ = -1 j 


(n) 
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= —3iy — 
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mi 

m2 

~4~~ 

2 


(3.25) 


which is the main result of this paper. For clarity pur¬ 
poses, we summarize the whole, new effective Hamilto¬ 
nian: 


Heff = 




(r^ + (Y) 


eff 


eff 


—fff + ffs^l (n X P) ■ Xi + ( —fff + ffs® j (nxp)-X2 

777-2 / \ 777i 


(0 ' O 


i-f 


1 


(i + ^) 


p^+f^-l] {n-pf 


^ ^2r2 - Af + + (n • x)ls^ ((«- x p) ■ xf^ej + Q4{p) 


1/2 


(3.26) 
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with 




(3.27) 


and with 


Af = 


94 41 

l.2u + 2uu^+l---n 


uu 




eff 


A®® = A®® (l + Qvu^ + 2(26 — ‘iv)vu^^ 


((n X p) ■ x)eff = 2(n • p){n ■ xo)(p • Ao) - P^(« • X)ls 

- (P • Xo) + (P^ - (« • P?) (X^)eff • 

(3.28) 

The A-potentials are 


where ^ 3 ^ denotes the (1,3)-Pade approximant taken with 
respect to the variable u = r“^. Notice that, as in 
Ref. [36], we do not take the Fade with respect to the 
variable r contained in 

We recall that the LO effective spin is defined as Xo = 
^Xi + ^X 2 , while the implementation of NLO spin- 
spin effects reads as follows: 


(x^)eff=X^ 




3 TOi \ 

2 7712 / 
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X 1 X 2 


(3.29a) 


{n-x)ls = in-Xof - 4 


21 9 mi\ n 9 

y 2 ^) ^ ^ 


21 9 7772 \ 9 9 


17 


77 


377 - 


9¥ 


P + TT - 


9 7771 

7 ^- 

4 7772 ! r 

9 7772 

7 ^- 

4 7771 j r 

37 / 7772 X ^ 

2 ¥ / r 


{n ■ Xif 

{n ■ X2f 
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(3.29b) 


The 25 coefficients that build the rather complex NLO 
spin-spin Hamiltonian in ADM coordinates get con¬ 
densed into 12 new contributions in the EOB. This re¬ 
sult is maybe not as striking as in the nonspinning sector 
(where, at 3PN order, the 11 ADM coefficients are re¬ 
duced to 3 EOB terms only), still it confirms the notable 
ability of the EOB in reproducing higher-order effects. 


IV. DISCUSSION: “FULL” AND “PARTIAL” 
INCLUSION 

One of the basic ideas behind the EOB is that of defin¬ 
ing a map between real and effective quantities. The 
complicated PN dynamics of spinning bodies, however, 
has forced a splitting of the concept of “effective spin” 
in the EOB. Instead of having one single map between 
{Si,S 2 ) and S^s, one has to distinguish between linear 
and squared spin, and even, as shown in this paper, be¬ 
tween different contractions of the spin vector with the 
dynamical variables, and map each component in a dif¬ 
ferent way. Having thus lost, so to say, the concept of an 
unique effective spin, one might ask up to which point 


r 


are we allowed to split the handling of it. The inclu¬ 
sion of NLO spin-spin effects we have presented in this 
paper obeys a simple rule: the spin mapping is only 
diversified when strictly necessary, in other words, the 
spin terms are as much as possible equally treated. We 
may call this the “full” approach, because of its formal 
and intuitive simplicity. One could also have followed 
another line of thought, pursuing the computational ef¬ 
ficiency rather than the conceptual unity. Such an al¬ 
ternative approach would consist in leaving untouched 
some spin terms of the effective Hamiltonian that actu¬ 
ally don’t contribute to the NLO spin-spin part of the PN 
expanded EOB. For example, the spin-orbit contribution 
Hso = AHso + P^Pi of the effective Hamiltonian only 
generates terms that are odd in the spins, and is thus 
unrelated to the spin-spin coupling. We may therefore 
apply the prescription ( |3.10 1 only to the orbital (even in 
the spins) part i7orb = a\/l -I- leaving Hgo as 

given by the “old” formulation (3.7a). We denote this 


kind of inclusion as “partial”. It is clear that this way of 
proceeding leads to an explicit expression for HeS which 
is shorter and simpler than in the “full” case, but it im¬ 
plies, at the same time, an additional and unnecessary 
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fragmentation of the effective spin. In addition, a prob¬ 
lem related to this approach is the fact that the “partial” 
Hamiltonian contains two different Delta potentials, Aj 
and A®®. Since these potentials are defined through a 
Fade approximant, which may contain poles, the “par¬ 
tial” Hamiltonian will in principle show twice as many 
poles than the “full” one, if the spin is nonzero. 

Let us finally remark that the “full” and “partial” 
Hamiltonians differ by terms that are odd (and at least 
cubic) in the spins, and therefore comparing them may 
help to distinguish the effects intrinsically related to the 
pure NLO spin-spin coupling from higher-order terms 
that get automatically resummed into the EOB after the 
NLO spin-spin inclusion. 


0.18 
0.16 
0.14 
O 0.12 

'3 0.10 
0.08 
0.06 




Figure 1. Reduced angular frequency (top) and binding en¬ 
ergy (bottom) at the LSO as a function of the effective spin 
parameter xo, in the case of equal masses and equal spins, for 
different EOB models. 

Fig. [2 shows the angular freguency Also and the bind¬ 
ing energy clso at the LSO (see also Refs. [51113Z]) as 
a function of the spin parameter xo for different EOB 
models, in the case of equal masses {ty = 1/4). The 
curves “3PN NLOssF”,“3PN NLOssP”, “4PN NLOssF” 
and “4PN NLOssP” denote the purely analytical EOB 
model discussed in this paper and in Ref. |36j . with 
NLO spin-spin coupling (“F” and “P” indicating the 
“full” or “partial” inclusion) and with next-to-next-to- 
leading order (NNLO) spin-orbit coupling (that is repro¬ 


duced by the gyro-gravitomagnetic factors calculated in 
Refs. [23134]). The label “3PN” or “4PN” refers to the 
purely orbital order, resummed with Fade or P 4 , re¬ 
spectively. The 4PN orbital order needs an additional 
term v {^al{iy) + (v) Inu) (Eq. (5) of Ref. [H]) in¬ 

side of At. The figure clearly shows that the effect of the 
4PN terms at the LSO is repulsive, since the frequency 
and the binding energy get increased. 

In Ref. [55J |57| , it has already been pointed out that a 
smaller effective spin squared leads to more bounded or¬ 
bits. It is thus not surprising that the “partial” inclusion 
is less bounded than the “full” one, as it can be seen in 
the figure. 

For completeness. Fig. [l] also shows two curves gen¬ 
erated by a different EOB model. Here, “tarl2” and 
“tarl4” denote the calibrated model of Refs. [551 155] . 
They both contain analytical information up to LO in the 
spin squared sector, with a calibration at the NLO level, 
and to NNLO in the spin-orbit sector, with calibration 
at next-to-next-to-next-to-leading order (NNNLO). The 
main difference between the two models is that “tarl4” 
reproduces the 4PN purely orbital order and is calibrated 
for a varying mass-ratio, while “tarl 2 ” only contains the 
3PN orbital order, and is calibrated for equal masses. 

Since the calibration is done at the waveform level, 
together with the tuning of some parameters related 
to the dissipative part, it is difficult to have a precise 
guess about the real accuracy of the curves “tarl 2 ” and 
“tarl4”. Hoping that the deviation from reality is not 
so large as to compromise a qualitative discussion of the 
plot, we can observe that the “partial” approach does 
not seem to show any particular advantage with respect 
to the “full” one. On the contrary, the “full” curves 
are generally closer to the corresponding calibrated ones. 
This might be interpreted as a further argument in favour 
of the “full” approach. Sistematically replacing, every¬ 
where in the EOB Hamiltonian, terms by their “effec¬ 
tive” equivalent is thus not only conceptually robust, but 
seems even to behave well numerically. 

A second interesting point is that the difference be¬ 
tween our EOB model and the calibrated “tar” mod¬ 
els is significantly smaller at the 4PN level than at the 
3PN. In particular, the frequency Also predicted by 
“tarl4” roughly follows the curves “4PN LSOssF”and 
“4PN LSOssP”, and lies in the gap between the two for 
Xo > 0.5. The maximal deviation between “4PN” and 
“tarl4” is 6.2% (“P”) and 8.7% (“F”), while between 
“3PN” and “tarl2” is 24.6% (“F”) and 26.7% (“P”). The 
same is true in the nonspinning regime. For xo = 0, the 
“4PN” value of Also deviates from “tarl4” by only 2.6%, 
whereas the difference between “tarl2” and “3PN” is of 
12.7%. 

In the case of the binding energy, there is a good corre¬ 
spondence for spin equal to zero (the deviations are 1.4% 
at 4PN and 5.0% at 3PN). However, for large spins, the 
difference is quite large (up to 39.2% (“F”) and 43.5% 
(“P”) at the 3PN level, and up to 20.9% (“F”) and 27.1% 
(“P”) at 4PN). Nevertheless, from 3PN to 4PN there is 
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still an improvement up to a factor ~ 2. This fact may 
strengthen the hope that the need for a calibration be¬ 
comes less urgent when higher-order analytical terms are 
included. 


V. CONCLUSION 

We have proposed a prescription for modifying the 
EOB Hamiltonian of Ref. [27], so that, when expanded 
in PN terms, it reproduces the correct NLO spin-spin 
coupling for general precessing orbits. This is a general¬ 
ization of the result of Ref. jSSj , where only equatorial or¬ 
bits with aligned spins had been taken into account. The 
implementation of the correct spin-spin terms is possible 
after a suitable canonical phase-space transformation of 
the ADM Hamiltonian. We have first shown the result 
in a rather general gauge, with 18 free gauge parameters 
(that actually reduce to 10 if we require symmetry un¬ 
der exchange of the particles). Then, a specific gauge is 
chosen. Being impossible, under the type of prescription 
we have considered, to impose a DJS-type gauge (i.e., to 
remove all new inclusions of type we have done 

the simple choice of looking for a gauge for which the 
maximal number of the coefficients entering the effective 
spin squared can be set to zero. It turned out that there 


is a unique gauge satisfying this criterion. 

In the end, a slightly different approach is taken into 
account, where the effective spin squared is left at LO in 
the spin-orbit sector of the effective Hamiltonian. This 
partial inclusion is less repulsive than the full one, and 
has the unpleasant feature of increasing the number of 
poles of the Bade approximant by a factor 2. A com¬ 
parison of the (gauge invariant) angular frequency and 
binding energy, taken at the LSO, with the calibrated 
EOB models of Refs. |35l [38], leads to the conclusion 
that the “partial” approach do not show any particular 
advantage with respect to the “full” one. By an Occam’s 
razor-like argument, the “full” approach, which is con¬ 
ceptually more consistent, may thus be preferable. As 
a last thing, the plot also brings the encouraging evi¬ 
dence that, as higher order terms are included, the purely 
analytical EOB approaches more and more a calibrated 
model even in the strong field. 
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